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Abstract. - We investigate the pseudo-Casimir force in a slab of material composed of 
nematically ordered long polymers. We write the total mesoscopic energy together with the 
jrt \ constraint connecting the local density and director fluctuations and evaluate the corresponding 

fluctuation free energy by standard methods, ft leads to a pseudo-Casimir force of a different 
type than in the case of standard, short molecule nematic. We investigate its separation 
'■^ , dependence and its magnitude and explicitly derive the relevant limiting cases. 

a 
o 

Q . 

Casimir effect is due to constrained fluctuations in media with long ranged correlations ||l| . 
The physical nature of the medium is not particularly important. Though the standard 
Q>^ ■ Casimir effect has been introduced for constrained electromagnetic field fluctuations [|| it 
CNj \ has been realized that other systems with long-range correlations exhibit a similar type of 
fluctuation driven interactions. Most notably critical fluids Q, smectic manifolds Q and 
Q^ I liquid crystals |5), all of them being prime examples of correlated fluids, give rise to a pseudo- 
^£. . Casimir effect which conies about through constrained thermal (as opposed to quantum) 
f ^ I fluctuations of order parameters. The nature of these order parameters of course depends on 
the system under study but they all exhibit massless fluctuation spectra that eventually lead 
C^ ' to long-ranged fluctuation interactions. 

Delimiting ourselves to the case of liquid crystals, the pseudo-Casimir effect has been re- 

' I searched in nematic, smectic and columnar uniformly ordered systems [^. In the case of 

^ • inhomogeneous or frustrated order as in the case of the hybrid-aligned cell characterized by 

Q I opposing surface fields or in the Freedericksz cell where frustration arises from competing bulk 

O ■ and surface fields Q , the interaction induced by director fluctuations is enhanced substantially 

by frustration, the enhancement being progressively stronger as the system approaches the 

, ■ transition from uniform to distorted structure. 

rN ' There seems to be some experimental evidence that backs up these theoretical predictions 0, 

j3 g] but there still lacks a definitive proof of the observability of the pseudo-Casimir interaction. 

■ ■ ■ It appears however that in the case of spinodal dewetting of 5CB on a silicon wafer |^ the 

pseudo-Casimir interaction is essential in giving a consistent interpretation of experimental 
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data, implying that this experiment can be regarded as the first observation of the pseudo- 
Casimir effect in Hquid crystals llfl] . 

In this contribution we shall investigate the pseudo-Casimir effect in the realm of confined 
polymers. As a first case study we will take a nematic polymeric liquid crystal such as one can 
observe in stiff polyelectrolytes above the isotropic - nematic transition, confined between two 
apposed planar surfaces. A physical realization of this case would be a nematically ordered 
DNA confined between two surfaces or simply cut to form a slab of macroscopically oriented 



sample 11 1. We believe that this case is particularly instructive since we will be able to 
connect different limiting results with the existing calculations for ordinary, i.e. short chain, 
nematics [pj. It will be shown that the polymeric nature of the nematogens gives rise to 
additional features of the pseudo-Casimir effect that distinguish it qualitatively from results 
derived in the case of short nematogens. Also we will show that in the case of nematic polymers 
the pseudo-Casimir interaction depends on the equation of state |ll| of the polymer nematic 
through its compressibility. 

The elastic deformation of an ordinary short molecule nematic |12|] with long range orien- 
tational order along the axis z, can be described with an average director n(r), with small 
fluctuations in the {x,y) directions: n(r) ~ {l,6nx:,Sny). In this case, if the splay, twist and 
bend elastic constants are the Frank constants Ki , K2 and K^ one obtains for the mcsoscopic 
elastic hamiltonian 

H„ = i / d'^rj^dz \Ki {V± ■ Snf + K^ (V^ x Snf + K3 {djnf] . (1) 

Besides fluctuations in nematic director one should also consider the fluctuations in local 
density of the molecules. For short nematogens the director and density fluctuations are 
decoupled and one need not consider this part of the mesoscopic free energy explicitly. For 
long, polymer nematics the situation is altogether different. In this case we have to consider 
the part of the free energy due to nonhomogeneous density fluctuations p(r) ~ po + Sp{r) in 
the Ornstein-Zernicke form 

Hp = i / d^^dz [b{Sp)^ + Be (V(5p)'] , (2) 

where B is the osmotic compressibility modulus and f is the density correlation length. Be- 
cause polymers are here considered to be infinitely long one also has to consider the fact that 
director inhomogeneities can relax only if accompanied by the simultaneous density relax- 
ation |l^. This leads to the constraint 

= V (pn) ~ djp + po (Vx • 5n) . (3) 

Adding all the components of the mesoscopic free energy and taking the constraint Eq. ^ into 
account via a Lagrange multiplier we obtain 

7^ = i / d\^dz \Ki {Vi_ ■ Snf + K2 (Vx x 6nf + X3 {d,5nf + 

+B{5pf + Be {V6pf + C (d.Sp + po (Vi • 5n))'] . (4) 

Writing 6n in the Helmholtz ansatz Sn — 6n^\ + Sn± we realize that the tangential part of 
the hamiltonian is decoupled from the constraint and its contribution to the fluctuational free 
energy should be the same as for short nematics 0. 

Obviously in the case of polymer nematics the longitudinal director and the density fleld 
fluctuations are coupled. The fluctuational free energy is now obtained via 

T = -kT \nZ = -kT \n V6n / VSp exp{-f3H), (5) 

J5n{d)=0 J5p{d)=0 
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where we have assumed that at the boundary of the sample the fluctuations in density and 
director field are quenched (i.e. 5n{d) — 0,6p{d) = 0). In principle one could also introduce 
the anchoring and tension energy of the bounding surfaces which would not bring any new 
qualitative features into our discussion but would make the computations less transparent. 

We now analyze the fiuctuational free energy of a slab of nematic material bounded at 
z = ±L/2 where the average director is aligned with the z axis. In the case of short nematogens 
(C = 0) the free energy decouples and we obtain two contributions. A long range pseudo- 
Casimir free energy due to director fluctuations and a short range term due to screened density 
fluctuations. Thus 



T{C = 0) 



kTS /K3 
l67 V ^ 



Ki 



C(3) kT ejL/O 



where S is the area of the surface, C(") tti6 Riemann's Zeta function and 

/•CJO 

e(-^/0 = / udu In (1 - e"") . 



(6) 



(7) 



Since L/^ ^ 1 and the function e(x) has an exponential behavior for large values of the 
argument the density fluctuations make a negligible contribution to the total free energy. Our 
result thus evidently reduces to the one obtained by Ajdari et al. Q. 

In the case that the director and the density fluctuations are coupled the calculation of the 
functional integral Eq. ra becomes more complicated. First of all we transform the hamiltonian 
Eq. by introducing 



Sp = poV± -a+Z 
The coupling term in this case reduces to 



Sn = —dzSi = —a. 



(8) 



C (d^p + PoV^ ■ Snf ^ Cpl [djf = C'/' 



and the hamiltonian is 



n^ ^Jd^r^dz 



(9) 



K3 2? + K, (Vxa)' + B (V^a + /)' + B^^ (V^lV^a + f)f 



Be{V^k+f] +C'.f 



(10) 



Clearly in the limit of strictly enforced density-director coupling C — > 00 and as a conse- 
quence / becomes independent of z, thus / = 0. In this case we can introduce the following 
linear transformation 

a — >a + ViA, \/\A^-f, (11) 

obtaining finally for the hamiltonian in the limit of strong coupling 



n^\Jd\^dz 



K^a^+Ki {V^-kf + B{V^-i^^ + Be iV^{V±-af + Be (V^-af 



(12) 



Since a is a longitudinal vector (Eq. m one obtains in the Fourier space 
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'^ = lT,Q Jdz 



i^saV 



K^ 



Be)Q%' 



5(1 



e'Q')Q'a^ 



(13) 



The mesoscopic hamiltonian obviously corresponds to a persistent oscillator, i.e. an os- 
cillator with an additional fourth order term in the derivatives. Since the different Fourier 
components are decoupled the evaluation of the functional integral proceeds straightforward 
and can in fact be reduced to a Feynman integral for second-derivative Lagrangian solved 
exactly by Kleinert ||lj.|l^]. Taking into account the boundary conditions for the fluctuating 
fields we obtain up to a multiplicative constant 



Z = 



{2-KkT)-'^ y/Ur[uJ^ \ujI - u:1\{ujI + t^i)^^ 



sinh(a;iiir3 I/)sinh(aj2-^3 L) 



^^^f^^icoshicoJ<-'/'L)cosHuj2K-'/'L) - 1)' 



where 



(14) 



(col 



= ^{Ki+Be)Q^K'i 



-1/3 



iL0iLU2f = B (1 + f g') Q^iiT; 



1/3 



(15) 



__ 1 /Q 

Introducing now ft± = K^ (cui ± UJ2) and with 



nliQ) 



Ki + Be 

2K. 



g±4AVi + e'Q' 



where 



(16) 



A = ^^ (17) 

Ki + Be ^ ' 

we obtain for the regularized fluctuation free energy per unit area, where the formally 
divergent bulk and surface free energy terms have already been subtracted 



7" = 



kTS 
An 



QdQ In 1 - 2--±e 



-Oi L 



cosh(J7„L) 






,-20+L 



(18) 



In the above equations fl^ can become negative for sensible values of the correlation length 
S,. In this case the coshr2_L has to be replaced by cos|ri_|L but the integral contains no 
dangerous divergencies and can be evaluated straightforwardly. By introducing 



T* 



kTS 

47r 



A^ and L* 



2{Ki + Be 



B 



(19) 



we realize that the dimensionless free energy T j T* obtained from Eq. Hq depends only 
on dimensionless separation LjL* and the dimensionless coupling constant 7 = A^. This 
coupling constant basically represents the ratio between the density correlation length ^ and 
the polymer nematic correlation length perpendicular to the average director A^^ p4|. 
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Fig. 1 - The dependence of the dimensionless free energy ^ on the dimensionless separation -^ 
obtained from Eq. hq for different values of the couphng constant 7. The two scaling regimes at small 
and large jr are clearly visible. 



In order to estimate the magnitude of the pseudo-Casimir forces in confined nematic poly- 
mers we now have to connect the macroscopic elastic constants Ki,K^, B and the correlation 
length ^ with the microscopic parameters of the systems. This is in general a very difficult 
undertaking. For DNA this type of analysis has been performed in ||ll|,|l^ and we just quote 
these results. 

If the interaction potential between the segments of the polymers is U{D), where D is 
the average spacing between the molecules perpendicular to the average director and the 
persistence length of the polymers is Cp, one has uW 



B 



U{D)ID 

kT Cp pi_ + U{D)/D 

V3 f d^ Id 



U{D)/Ct 



(20) 



One should realize here that since DNA is a polyelectrolyte the interaction potential U{D) 
depends also on the ionic strength of the bathing solution. 



This scaling of elastic constants leads to the conclusion that in Eqs. |2l| and 22 we have 
■j^ "^ i and ^ ^ 1. Thus the larger the D, i.e. the smaller the density of the polymers, 
the larger this ratio becomes. The magnitude of the pseudo-Casimir interaction thus de- 
pends fundamentally on the equation of state of the nematic polymers through the osmotic 
compressibility B and the density perpendicular to the average nematic director, p±. 

Fig. 1 presents a plot of dimensionless free energy as a function of dimensionless separation. 
The curves are plotted for various values of the coupling constant 7 which is plotted as a func- 
tion of the average spacing between polymers on Fig. 2 for the case of DNA. One can clearly 
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Fig. 2 - The dependence of the magnitude of the pseudo-Casimir interaction T* , the characteristic 
length L* and the coupling constant 7 on the average separation between the polymers obtained 
from |pl| in the case of 0.5 M bathing ionic solution. Note that 7 and J-'* have the same scale. L* 
saturates for small values of DNA density (large D) at about 140 A. 



discern two regions with approximately L~^ and L~^ behavior separated by a characteristic 
spacing Lq — C(j)L* , where C(7) is a numerical factor depending on the couphng constant 
that can be read off Fig. 1. The numerical evidence shows that the caracteristic spacing Lq 
grows approximately linearly with the average separation between the polymers D, ranging 
from about 15 to 20 nm aX D — 2 nm to around 150 nm at D = 8 nm, which is between one 
and two correlation lengths S,- 

Fig. 2 shows the dependence of the magnitude of the pseudo-Casimir interaction T* , the 
characteristic length L* and the coupling constant 7 on the average separation between the 
polymers obtained from the measured equation of state for DNA in the case of 0.5M bathing 
electrolyte solution (for details see ||ll| ). 

One would of course like to have also approximate analytical formulae for the pseudo- 
Casimir interaction Eq. |l^. There are only two limiting cases where the integral in Eq. |8| 
can be evaluated analytically, the large separation and the small separation regime. First of 
all for small L the dominant behavior of the integral is obtained at large Q leading in this 
case to 

HL « io) = - ^ ( ^, -^ + ^ 1 C(3). (21) 



16ttL^ 



/fi(l + 4Ae)- 1 + ^e 



Ki 



The other analytically tractable limit is in the case of large L where the dominant contri- 
bution to the integral comes from small Q behavior of the integrand 



T{L » Lo) = - 



kTS 
TGtT 



2X3C(5) , K^m 



BL'i 



KiL^ 



(22) 



J. DoBNiKAR and R. Podgornik : Pseudo-Casimir interactions 



where the length scale Lq can be read off the graph on Fig. 1 for various values of 7 (see above). 
Thus for small L the confined polymer nematic behaves basically as a standard nematic M 
with a renormalized magnitude of the fluctuation interaction. For large L the dependence 
is altogether different, deviating essentially from the standard nematic case. The crossover 
between the two regimes moves to larger L as the polymers become less dense (see Fig. 2) 
but eventually saturates for small polymer densities. 

The pseudo-Casimir interaction for long nematic polymers thus differs qualitatively from 
the one obtained in the case of short nematics M. Though it still decays algebraically with 
L, and thus qualifles as a long range force, it decays faster (i"**) then in the case of short 
nematics {L~'^). Nevertheless its absolute magnitude, depending on the equation of state 
of the nematic polymer and the regime of spacings L, can become comparable to and even 
larger than in the case of short nematics with the same Frank constants. The L""* behavior 
stems essentially from the polymer stiffness and should thus be a salient feature of the pseudo- 
Casimir interactions whenever the mesoscopic hamiltonian contains the squares of higher then 
the first derivatives in the order parameter. These type of systems will be studied in our future 
work. 
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